Concentration-dependent mobility in organic field-effect transistors probed by 
infrared spectromicroscopy of the charge density profile 
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We show that infrared imaging of the charge density profile in organic field-effect transistors 
(FETs) can probe transport characteristics which are difficult to access by conventional contact- 
based measurements. Specifically, we carry out experiments and modeling of infrared spectromi- 
croscopy of poly(3-hexylthiophene) (P3HT) FETs in which charge injection is affected by a relatively 
low resistance of the gate insulators. We conclude that the mobility of P3HT has a power-law den- 
sity dependence, which is consistent with the activated transport in disorder-induced tails of the 
density of states. 
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Transport characteristics of organic semiconductors 
are important for optimizing the performance of organic 
FETs. Due to substantial disorder, polaronic, and in- 
teraction effects, transport mechanisms in organics are 
complicated and standard model assumptions used for 
inorganic devices may be unrealistic. For example, the 
canonical model of the FET [lj assumes that the mobil- 
ity \i is independent of the areal carrier density N in the 
accumulation (or inversion) layer, while theoretical and 
experimental investigations of organic conductors point 
to the contrary [2|, y, U, [5| . Unfortunately, the extraction 
of mobility from traditional transport measurements is 
complicated by the nonlincarity of the contact resistance, 
which may obscure the intrinsic mobility of the organic 
semiconductor [1, 0|- 

Recently, we have reported a direct imaging of the 
charge distribution in the accumulation layer of organic 
FETs using infrared (1R) spectromicroscopy |8|. In this 
Letter we show that this contactless technique is com- 
plementary to traditional DC transport measurements. 
Specifically, the carrier density profile in the conducting 
channel obtained from the IR scans can be used to extract 
the intrinsic behavior of the density-dependent mobility, 
avoiding artifacts due to the contacts. 

As a case study, we investigated P3HT-based FETs 
in which an inhomogeneous charge density profile was 
formed due to a non-negligible leakage current through 
the gate dielectrics. Below we show that the idealized 
model [1] fails to correctly account for the charge density 
distribution observed in these FETs. Instead, we find a 
good agreement with the model of a density-dependent 
mobility, 
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FIG. 1: (a) FET schematics (b) A sketch of the V-shaped 
electrodes between which the IR scans were taken, (c) The 
voltage-induced absorption spectrum Sad of one of the FET 
devices studied. The sharp absorption peaks are the IR-active 
vibrational modes of P3HT. 



which is derived theoretically below assuming transport 



is activated and the density of states has an exponential 
tail with the characteristic energy width E*. Since E* 
depends on disorder, exponent (3 is not known precisely. 
From the fits to our IR spectromicroscopy data we find 
1 < (3 < 4 at temperature T — 300 K, where all the 
measurements are done. This implies E* ~ 0.1 eV, which 
agrees with other experiments on P3HT Q. The power- 
law behavior |T]) has been previously observed both in 
P3HT [5] and in other organic conductors {2J. Thus, it is 
common to this class of materials. 

Let us proceed to the description of our experiment. 
The structure of the FETs is shown in Fig. QJa). In 
these bottom-contacted devices, source and drain Au 
electrodes were deposited on Si02/Ti02/n-Si substrates 
followed by the spin coating of a 4-6 nm thick P3HT film, 
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which served as the electronically active material. To 
examine length scales associated with charge injection, 
we fabricated devices with a V-shaped electrode pattern 
[Fig. Gib)] where the electrode separation varied from 
0.01mm at one side to 4 mm at the other. To image 
the injected charges we have employed the IR spectro- 
microscopy, which involved scanning a focused IR beam 
50-100 pm in diameter across the region between the 
V-shaped electrodes and monitoring the changes in the 
P3HT absorption spectrum at a given gate voltage [H . 

The spectroscopic "fingerprints" of the injected 
charges in P3HT are the IR-active vibrational modes 
[Fig. QJc)] and the polaron band (a broad absorption 
maximum at u ~ 3500 cm -1 , cf. Fig. 2 of Ref. 0]). From 
the oscillator strength of the former, we computed the 
gate-induced charge density N(x) as a function of the dis- 
tance x from the source electrode. (These results agreed 
within the experimental accuracy with iVYx) computed 
from the polaron band spectra, see Ref. [8| for details.) 
In Fig. [21 we plot N(x) normalized to its value at x = 
for three representative devices. Our main result is that 
the injected charge density decays with the distance from 
the source and vanishes at x ~ 800 /im. Below we show 
that the spatial variation of the charge density contains 
valuable information on the iV-dependence of the mobil- 
ity in accumulation layers of organic FET. 

Indeed, in the experiment a non-negligible leakage cur- 
rent ji through the dielectric [lj| was detected. The 
corresponding effective resistivity of the dielectric layer 
pi = V/ji was a function of V, the voltage difference be- 
tween the polymer and the gate. However, pi « const 
was observed [l(J at V > 5 V and, for simplicity, we 
assume this to be the case. In the gradual-channel ap- 
proximation suitable for our FETs we obtain the balance 
of surface (j s ) and leakage (ji) currents as follows: 



Vj s = V (ctW) = V/ pi 



(2) 



For the sheet conductivity a we adopt the form a — 
&a + 0V, where ay and ay are the sheet conductivities of 
the accumulation layer and the residual mobile charges, 
respectively. Initially, let us make the conventional as- 
sumptions that (i) fx = cr a /N = const (ii) ay <§; ay, and 
(iii) V = eN/C, where C = const is the capacitance per 
unit area (C ~ 0.14 /iF/cm 2 in our FETs). This entails 



V(iVViV) ~ N/ienpi) 



(3) 



Since the characteristic length of charge inhomogeneity 
is smaller than the lateral extension of the source and 
drain electrodes in the V-shaped region [Fig. [2(b)] , it is 
permissible to treat them as infinite metallic half-planes, 
x < and x > L, respectively. In this case, all vari- 
ables depend only on x G (0. L). Suppose first that L is 
very large (the source and drain are far apart), then the 
physically relevant solution of Eq. ([3J) is 



Interestingly, it predicts that the charge injection has 
a definite cutoff length ay. This occurs because of the 
quadratic nonlinearity of Eq. ([3]) . Treating ad- 
justable parameter, we plot Eq. (0| in Fig. [2] (the curves 
labeled (3 — 0). Alternatively, we can estimate ay based 
on a and pi deduced from the DC transport measure- 
ments (HE!, cr - lO- 7 ^- 1 and pi ~ 10 5 ftcm 2 , which 
yields ay ~ 2500 pm. Comparing with Fig. [21 we see that 
the observed ay is significantly smaller. Furthermore, in- 
stead of a more gradual decrease in the conventional the- 
ory, in experimental data N(x) exhibits an abrupt drop. 
This suggests that the degree of nonlinearity of the actual 
current balance equation ([2]) is higher than quadratic. 

The field-dependence of p is probably not the main 
source of this effect because we study P3HT of a rather 
high mobility p = 0.05-0.12 cm 2 /Vs, where non-Ohmic 
behavior is not as pronounced as in low-mobility organ- 
ics 0,0,0,01) an d also because T was relatively high 
On the other hand, there is a natural reason for the iV- 
dependence of p [conveyed by Eq. (JTJ]. At densities we 
deal with, N < 10 13 cm~ 2 , the system approaches the 
insulator-to-metal transition region [I3 |. The chemical 
potential £ presumably resides in the density-of-states 
tail of disorder-induced traps but at temperatures of in- 
terest such traps are shallow (a few ksT deep). There- 
fore, transport is dominated by the activation of carriers 
to either the mobility edge or the so-called transport en- 
ergy level [11,111 and it is characterized by a modestly large 
activation energy i.e., ay oc exp^/Zc^T) (the mobil- 
ity edge is taken to be the energy reference point) [HI]. 
Because of this exponential relation between the conduc- 
tivity and C, the dependence of £ on N is amplified in 
a a , producing the sought nonlinear effect. 

It is possible to express ay directly as a function of 
N if the density of states g(E) is known. In conducting 
polymers g(E)is often modeled by Gaussian or expo- 
nential forms [3j, which are observed experimentally [9|. 
In either case we obtain N(£) = f dEg(E)f(E,() ^ 
N t exp(C/-E*), up to logarithmic factors. Here f(E,Q 
is the Fermi-Dirac distribution function, E* > ksT is 
the characteristic energy width of the band tail, and N t 
is the total number of states per unit area in this tail. We 
see that C - ln(N/N t ), leading to ay ~ <r*(N/N t y +1 
with (3 + 1 = E^/ksT. Consequently, the mobility is 
p(N) = (T a /eN cx N 13 , in agreement with Eq. |T]). We 
note that a nearly identical p(N) dependence also follows 
from the hopping transport model [! i]. Equation (JJ) 
was also postulated in Ref. 0- Substituting a cx N l3+1 
into Eq. ([2]) but still neglecting ay, we obtain 



N ~ N(0) 1 



.-i+i 



y/2 g + 3) a(0) Pl 
13 + 1 



(•5) 



N(x) ~JV(0)[1- (x/ay)] 2 , x < ay = y/6a{0) Pl . (4) 



with Eq. ([4]) being recovered for j3 = 0. Three repre- 
sentative curves computed according to Eq. ([5]) are plot- 
ted in each panel of Fig. [21 Clearly, for all of the three 
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yields the algebraic equation for N(x), 
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FIG. 2: Charge density profiles measured by the IR spectro- 
microscopy in three nominally similar devices. The solid lines 
are theoretical fits (see text). 



FETs studied, one can find /3s somewhere in the range 
I < f3 < 4 which give a significantly better fit to exper- 
iment than (3 = 0, both in the functional form and in 
the value of a;* . Apparently, Eq. ([T]) is more suitable for 
modeling the charge injection profile in a "leaky" FET. 

For future reference we mention that a refinement of 
the model including the residual conductivity, 



a*(N/N t )^ 



ay. 



(6) 



predicts the injected charge density profile of the form 



— = -V^T^ + tt^J ln /—r- r + const ' ( ? ) 
a;* p + A y/u + a — y/a 

u = [N(x)/N(Q)f+\ a=(fi + 3)<7 r /[2o- a (0)]. (8) 

Accordingly, N(x) oc exp(— x/y/cr r pi) at x > x*, i.e., it 
has the exponential decay instead of the abrupt threshold 
predicted by Eq. ([5]). However, since a r /a(0) ~ 10~ 4 - 
10~ 3 (the on-off ratio) the pre-exponential factor of 
this tail is very small, producing apparently no signature 
in our experiment, see Fig.O 

Finally, let us also discuss the L < x* case where the 
injected charge is able to reach the drain (as in SiC>2 
devices we studied previously Here N(x) depends on 
the bias Vsd = V(0) — V(L). Also, the right-hand side 
of Eq. ([2]) can be set to zero because we can normally 
neglect the leakage current ji compared to the source- 
drain current j s . Equation ^ remains solvable 1J| and 



a*Nt (N_ 
(3 + 2\N t 



(3+2 



a r N = — — (x„ — x), (9) 
e 



where x p is the integration constant. At V(L) — > 0, 
we have x p — > L and, except very near the drain, the 
first term on the left-hand side dominates. In this case 
the charge profile is very similar to that of the injection 
case, Eq. (|S|), except a;* is replaced by L and the power- 
law exponent 2/(/3 + 2) is changed to a smaller number 
l/(/3 + 2). Although we were not able to verify this for- 
mula by our technique (in part, due to a finite spatial 
resolution ~ 100 /im) we note that Biirgi et al. have 
studied this regime in similar FET devices using scanning 
potentiometry. These authors were able to fit their data 
to the formula obtained by setting ay = and (3 = 1 in 
our Eq. © 0. 

In conclusion, when the charge profile in the or- 
ganic semiconductor is strongly inhomogeneous, as in the 
charge injection [8j or the current saturation [5( regimes, 
its modeling should include the density dependence of the 
carrier mobility. Our experiment is consistent with the 
power law ([T]), which has a natural physical motivation. 
Our measurements and the method of analysis are free 
from problems introduced by the contact resistance. 
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